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[0001] The field of the invention relates to medical devices, and more particularly, to methods 
for tracking catheters such as those used to conduct ultrasonic imaging. 
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Background of the Invention 



5 [0002] In the medical arts, catheters are frequently used to diagnose and treat various disorders 
in a patient, such as clogged or blocked blood vessels. A catheter is introduced into a blood 
vessel of a patient by, for example, making an incision in the patient over the blood vessel and 
inserting the catheter into the blood vessel of the patient. A catheter operator such as a physician 
then maneuvers the catheter through the blood vessels of the patient until the catheter is properly 
10 situated to diagnose or treat the disorder. Similar techniques are used to insert catheters into 
other types of lumens within a patient. 

[0003] In maneuvering the catheter through the blood vessels or other lumens within the patient, 
there is a recurrent need to know the location of the catheter within the body space of the patient. 
Conventional imaging systems create an image of the blood vessel or other lumen which make 

15 the lumen appear as a straight tube, and provide no concept of 3 -dimensional ("3-D") spatial 

relationships. In the patient, however, the lumens curve about, and contain branches that branch 
off at various angles from the lumen. If the position in three dimensions of the imaging head on 
the catheter can be determined, then through use of three-dimensional imaging software, the true 
positions and locations of the curves, twists, and turns, as well as the locations of the branch 

20 points, of the lumens can be determined. Knowing the true positions allows a more accurate map 
of the patient to be created, which yields more effective diagnosis and treatment of the patient. 
For example, gathering accurate 3-D position data allows for an accurate blood flow map and 
consequent blood flow monitoring and modeling. 

[0004] Traditionally, X-ray technology has been used to provide a global roadmap of X-ray 
25 visible devices, showing their position within the patient. However, an X-ray image, being a two- 
dimensional projection, can only provide partial information on the 3-D shape of the catheter 
path. Furthermore, prolonged exposure to X-rays may be harmful to the patient, and it is 
therefore desirable to avoid such exposures. Thus there is a need for a tracking system which 
can easily determine the location of a catheter within a patient, without exposing the patient to 
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harmful side effects, and which can be used with a wide variety of catheters or other imaging 
medical devices. 

[0005] To overcome the problems inherent with X-ray tracking of catheters, various 
technologies have arisen which attempt to gather positional information about the location of a 
5 catheter within the patient, without the harmful side-effects of X-ray technology. Among such 
technologies are tracking systems which gather positional information using electromagnetic, 
optical, mechanical, acoustic, and/or inertial sensing elements. Many of these technologies 
require the addition of extra elements to the catheter, to allow it to be tracked within the patient. 

[0006] Therefore there is a need for an improved method of tracking catheters. 

10 Summary of the Invention 

[0007] For imaging catheters, the disadvantages of X-ray tracking can be avoided, without 
needing any additional equipment added on to the catheter, by relying on the data contained in 
the images collected by the imaging catheter itself to determine the position of the catheter 
within the body. This improved method can also be used with other forms of catheters, as long 
15 as the catheter has some ability to gather data about its immediate surroundings. 

[0008] In this method, a first image gathered by the imaging catheter is compared to a second 
image gathered by the imaging catheter, and this comparison is used to compute one or more 
offset angles between the first and second images. This data is used to determine the relative 
position of the second image with respect to the first image. By making these determinations for 

20 each of a series of images, the orientation of the entire series of images, in three dimensions, can 
be determined. Since the orientation of an imaging catheter image is determined by the 
orientation of the imaging element at the tip of the imaging catheter, this method allows the 
position of the imaging element to be determined. This method also allows an imaging system to 
create a true or more true three-dimensional representation of the lumen that the catheter is 

25 traveling through. 

[0009] Other systems, methods, features and advantages of the invention will be or will become 
apparent to one with skill in the art upon examination of the following figures and detailed 
description. It is intended that all such additional systems, methods, features and advantages be 
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included within this description, be within the scope of the invention, and be protected by the 
accompanying claims. 

Brief Description of the Drawings 

[0010] In order to better appreciate how the above-recited and other advantages and objects of 
5 the present inventions are obtained, a more particular description of the invention briefly 
described above will be rendered by reference to specific embodiments thereof, which are 
illustrated in the accompanying drawings. It should be noted that the components in the figures 
are not necessarily to scale, emphasis instead being placed upon illustrating the principles of the 
invention. Moreover, in the figures, like reference numerals designate corresponding parts 
10 throughout the different views. However, like parts do not always have like reference numerals. 
Moreover, all illustrations are intended to convey concepts, where relative sizes, shapes and 
other detailed attributes may be illustrated schematically rather than literally or precisely. 

[0011] FIG. 1 is an imaging system in accordance with an example embodiment of the invention. 

[0012] FIG. 2 is a series of image slices of a volume of material. 

15 [0013] FIG. 3 is a series of image slices of a volume of material, including example objects 
within the volume of material. 

[0014] FIG. 4 is the series of image slices, showing the cross-sections of the example objects 
within each image slice. 

[0015] FIG. 5 is an example of a graph of the loss of correlation between two pairs of planes. 

20 [0016] FIG. 6 shows the fitted exponential functions for the graph of FIG. 5. 

[0017] FIG. 7 is a graph of the relationship between derivatives of exponential functions and 
angles of separation. 

[0018] FIG. 8 shows the position of an image plane in polar coordinates, with respect to a 
reference plane. 

25 [0019] FIG. 9 shows the position of a second image plane in polar coordinates, with respect to a 
reference plane. 

Detailed Description of the Preferred Embodiments 
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[0020] Turning to FIG. 1, an imaging system for use with an example embodiment of the 
invention includes an imaging catheter 10 coupled to a computer 20 which is executing imaging 
software to display images on a video display 30. The imaging catheter 10 includes an imaging 
head 12 located at the tip of the imaging catheter 10. The imaging head 12 may be an ultrasonic 
5 transducer, a light based imaging device, or any other device which gathers imaging data. The 
imaging catheter 10 is inserted into a lumen in a patient, such as a blood vessel. The imaging 
head 12 gathers image data from within the lumen of the patient, for example by collecting a 
cross section image of the lumen. In one specific embodiment, the imaging head 12 rotates 
about the axis of the catheter 10, and gathers a circular band of image data from the lumen by 

10 reflecting sound waves off of the lumen. Alternatively, light waves can be reflected off of the 
lumen and received by the imaging head 12. The imaging head 12 sends this image data down 
the catheter 10 along a data communication channel 15 such as a wire or fiber optic channel, to 
the computer 20. The computer 20 receives the image data, optionally processes the image data, 
and then displays the image data on the video display 30. In alternate embodiments, other image 

15 display devices may be used in place of the video display 30, such as holographic display units, 
printers, projectors, etc. In this example embodiment, the images taken by the imaging head 12 
represent cross-sections with respect to the axis of the lumen in which the catheter 10 is 
traveling. In alternate embodiments, images captured by the imaging head 12 represent slices 
taken at angles other than cross-sections with respect to the axis. 

20 [00211 Since the lumen is twisting, bending and curving about inside the patient, the lumen axis 
is constantly changing orientation, and thus the orientation of the images is also constantly 
changing. Additionally, the catheter operator may elect to alter the orientation of the imaging 
head within the lumen, for example to take an image of a particular portion of the lumen at a 
different angle in order to gather additional data about the lumen. To provide an accurate 

25 determination of the position of the imaging head 12, and thus provide an accurate map of the 
lumen, it is useful to determine the relative offset of each image from the image captured just 
previously. To conserve computing resources and time, it is possible to skip some of the images, 
and compute the relative offset of an image with some image other than the image captured just 
previously. This may result in a corresponding decrease in accuracy, but this may be acceptable 

30 depending on the particular situations the catheter 10 is being used in. 
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[0022] To simplify the discussion and more clearly explain the method of an embodiment of the 
invention, the following disclosure of the process of image comparison and position 
determination will use the example image slices shown in FIG. 2. The series of intersecting 
planes lOOa-j within the volume box 110 represent image slices taken of a volume image of 
5 some biological material, such as a segment of a blood vessel or other lumen within a patient. 
The intersecting planes lOOa-j all meet at the line of intersection 105, down the left side of FIG. 
2. The intersecting planes lOOa-j could be, for example, a series of images captured by the 
imaging head 12 as the catheter 10 was guided around a bend in the lumen, beginning with the 
plane 100a at the bottom of the volume box 110, and progressing through the other planes lOOb-j 
10 towards the top of the volume box 110. 

[0023] Each plane lOOa-j contains a slice of image data, such as an ultrasound image, or a light 
based image. The image data contained in each plane lOOa-j changes, as the planes progress 
through the volume box 110. For example, turning to FIG. 3, the volume box 1 10 contains a first 
shape 120 and a second shape 130. Actual ultrasound or light based images would be more 

15 complex than the example presented here, but the same principles as disclosed apply to any sort 
of image data. This example has been simplified for clarity. The first shape 120 begins at the 
bottom of the volume box 1 10, at plane 100a, and continues up through each of the planes 100b- 
j, ending in the space above the last plane lOOj. The second shape 130 is the same size and 
orientation as the first shape 120, but in a different location in the volume box 110. The second 

20 shape 130 begins at the intersection point with the plane 100b, and extends upwards through 
each of the planes lOOc-f, ending in the space between plane lOOf and plane lOOg. Planes 100a 
and lOOg-j do not intersect the second shape 130. 

[0024] The planes lOOa-j intersect the shapes 120, 130, capturing a slice of image data for each 
shape 120, 130 at the intersection. The resulting series of planes lOOa-j containing the slices 

25 120a-j, 130a-e of the shapes 120, 130 are shown in FIG. 4. Plane 100a contains a small slice 
120a of the first shape 120. Plane 100b contains a slightly larger slice 120b of the first shape 
120, and a small slice 130a of the second shape 130. The slices 120a-j of the first shape 120 
become progressively wider as the planes progress through the slices 120a-j of the first shape 
120 from plane 100a to plane lOOj. The slices 130a-e of the second shape 130 become at first 

30 wider and then narrow back down and eventually cease, as the planes progress through the slices 
130a-e of the second shape 130. 
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[0025] As can be seen in FIG. 4, there are differences in the image slices on the planes lOOa-j, as 
the planes progress through the shapes 120, 130. However, these differences are more 
pronounced on the right side of each plane lOOa-j (for slices 130a-e of shape 130), as compared 
with the differences on the left side of each plane lOOa-j (for the slices 120a-j of the same sized 
5 shape 120). These differences occur because the separation distance between adjacent planes, 
such as planes lOOa-b, is greater at the right side than it is at the left side. Since the two planes 
lOOa-b are further apart from each other on the right side, the correlation between the image 
slices 130a-b on the right side is correspondingly lower than the correlation between the image 
slices 120a-b in the same two planes lOOa-b on the left side. 

10 [0026] In the example of FIG. 4, the differences are determined by making a high-level 

comparison of the appearance of the two objects 120, 130 on each plane being examined. For 
other types of data, other ways of determining the difference and the variation of the difference 
across the image are possible. For example, if the data represents the density of a particular 
point on the plane (such as an ultrasound image), then the difference is calculated by measuring 

1 5 the difference in density values for the same point on each of the two planes. If the data 

represents visual image data, then the difference may be calculated by measuring the difference 
in some value or property associated with the visual image data, such as the color, hue, 
saturation, reflectivity, or density of the image at the same point on each plane. 

[0027] In general, the correlation between any two of the planes lOOa-j is greatest at the line of 
20 intersection 105 along the left side of each plane lOOa-j, where all of the planes contain exactly 
the same image data. Correlation between pairs of planes lOOa-j is gradually lost as we progress 
across the planes lOOa-j from left to right. We will use the term "correlation loss" or "loss of 
correlation" to describe the local difference between the images being compared, which may be 
computed by a variety of methods. The particular method of computing the difference between 
25 images is a design choice for those skilled in the art, and is not critical to the embodiments of the 
invention disclosed herein. In this example, for ease of explanation, the correlation loss was 
measured across the planes lOOa-j from left to right. The correlation loss may be measured in 
any direction across the planes lOOa-j without affecting the results of the methods discussed 
herein. 
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[0028] An example considering images taken across two pairs of planes (all intersecting along a 
common line) demonstrating the behavior of image differences with increasing distance from the 
common intersecting line is shown in the graph in FIG. 5. The horizontal axis 140 of the graph 
of FIG. 5 measures the distance away from the intersection line 105. The vertical axis 150 is a 
5 measure of the difference between the two planes lOOa-j being compared. The lower data line 
160 represents sample data gathered by comparing the image data on plane 100a with the image 
data on plane 100b. The upper data line 170 represents sample data gathered by comparing the 
image data on plane 100a with the image data on plane 100c. It is evident that the images in the 
pairs being considered start out being highly correlated (as evidenced by the small difference 
10 between them) and then gradually lose correlation (as evidenced by the difference becoming 
larger). The rate of correlation loss between planes 100a and 100b is less than the rate of 
correlation loss between planes 100a and 100c. This is due to the fact that the angle between 
plane 100a and plane 100b is smaller than the angle between plane 100a and plane 100c. 

[0029] Turning to FIG. 6, the difference between any two planes lOOa-j can be approximated by 
15 an exponential function of the form f(x)=A(l-e _Xx ). Here A represents the value ultimately 

attained by the function and X is a parameter describing the rate at which the function increases 
as it attains the value A. The particular exponential function representing the behavior of the 
difference between any two planes lOOa-j may be approximated using a best fit approximation 
method. The lower fitted exponential function 165 is an approximation for the lower data line 
20 160. The upper fitted exponential function 175 is an approximation for the upper data line 170. 

[0030] Calculating the derivative of each fitted exponential function 165, 175 at the origin (0,0), 
yields the value AX for each fitted exponential function 165, 175, which value is a good 
approximation of the angle of separation, in the direction of the correlation loss comparison, 
between the two planes lOOa-lOOj being compared. Thus, the derivative of the lower fitted 
25 exponential function 165 at the origin is an approximation of the angle of separation between 
planes 100a and 100b. Similarly, the derivative of the upper fitted exponential function 175 at 
the origin is an approximation of the angle of separation between planes 100a and 100c. 

[0031] FIG. 7 shows a graph of the relationship between 1) the derivatives of a series of fitted 
exponential functions approximating sample correlation data comparing each of planes lOOb-j 
30 with plane 100a and 2) the angle of separation 8 between each of planes lOOb-j and plane 100a. 
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For this sample data, the angle of separation 0 between each pair of planes was already known. 
The horizontal axis represents the AX values which are approximated from the correlation 
comparisons discussed above. The vertical axis represents the known angle of separation 0. The 
series of data points 180a-j represent the relationship between the approximated AX values and 
5 the known angles of separation 0 for each plane lOOa-j as compared with the reference plane 
100a. Data point 180a represents the angle of separation between plane 100a and itself (i.e. 
zero), as related to the approximated AX value for the correlation comparison of reference plane 
100a with itself (i.e. zero). Data point 1 80b represents the angle of separation between reference 
plane 100a and plane 100b, as related to the approximated AX value for the correlation 

10 comparison of plane reference 100a with 100b. Similarly, the data points 180c-j represent the 
angle of separation between the respective planes lOOc-j and the reference plane 100a, as related 
to the approximated AX value for the correlation comparisons of the respective planes lOOc-j 
with the reference plane 100a. In this example, the fitted line 185 demonstrates that the 
relationship between the AX values and the angles of separation 0 is a linear relationship, and that 

15 the AX values are proportional to the angles of separation 0. Thus the proportionality constant 
can be computed by dividing the known angle of separation values by the AX values. Therefore, 
using the process just described the angle of separation between any two planes in any given 
direction can be estimated by analyzing the rate of loss of correlation between the two planes in 
that direction. 

20 [0032] Turning to FIG. 8, once the angle of separation 0 between any plane lOOb-j and the 
reference plane 100a is known, principles of geometry allow the relative position of the plane 
lOOb-j with respect to the reference plane 100a to be calculated. For example, referring to plane 
lOOg, the plane lOOg is defined by the two lines 105, 190. Since both of these lines are 
calculatable with respect to the reference plane 100a using data gathered from the two images on 

25 the planes 100a and lOOg, the position of the plane lOOg is calculatable with respect to the plane 
100a, using data gathered from the two images on the planes 100a and lOOg. 

[0033] Assume that the reference plane 100a is defined to begin at (0,0,0) in a three-dimensional 
polar coordinate system (p, 0, <£), with the reference plane 100a defined by the two intersecting 
lines (0,0,0)- (p, 0, 0) and (0,0,0)- (p, 0, tt/2), where p represents the width of the image on the 
30 reference plane 100a, 0 represents the angle of inclination above the reference plane 100a, and <j> 
represents the angle of counterclockwise rotation from the origin line (0,0,0)- (p, 0, 0). The 
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position of the plane lOOg is defined by the position of the intersecting line 105, represented in 
polar coordinates as (0,0 5 0)- (p, 0, 0) and line 190, represented in polar coordinates as (0,0,0)- (p, 
0, 0). The value for p is known from the dimensions of the reference plane 100a, and the value 
for 8 is calculated using the correlation comparison of the image data in the plane lOOg as 
5 compared with the reference plane 100a as discussed above. Therefore, the position of the plane 
lOOg is determined using the information known about the reference plane 100a and the 
information gathered from the image data correlation comparison 

[0034] Note that the example presented above limited the correlation comparison to planes 100a- 
j that only varied from each other in one direction, since all of the planes lOOa-j shared the line 

10 of intersection 105. This was done for ease of explanation of the principles of operation of the 
improved method of tracking a catheter using image data, but alternate embodiments of the 
method can easily calculate the position of other planes which have any arbitrary relationship to 
a reference plane 100a. For example, turning to FIG. 9, a plane 100k varies from the reference 
plane 100a in two directions. There is an angle 9\ between the reference plane 100a and the 

15 plane 100k, where the two planes 100a, 100k intersect the (p, 0, 0) axis, and an angle 6 2 between 
the reference plane 100a and the plane 100k, where the two planes 100a, 100k intersect the (p, 0, 
7r/2) axis. 

[0035] To determine the angle 0\ between the reference plane 100a and the plane 100k, the 
correlation loss rate in the direction (0,0,0) to (p, 0, 0) is computed and approximated to an 

20 exponential function as described above. To determine the angle 62 between the reference plane 
100a and the plane 100k, the correlation loss rate in the direction (0,0,0) to (p, 0, tc/2) is 
computed and approximated to an exponential function as described above. Note that while the 
correlation loss began at zero for the example of FIGS 1-8, this will not always be the case for 
the example of FIG. 9. For any sampling line (such as the sampling lines 200, 210, 220, 230) 

25 across the plane 100k other than the sampling lines following the axes from (0,0,0) to (p, 0, 0) 
and from (0,0,0) to (p, 0, 7r/2), the correlation loss function will begin with a non-zero value, 
indicating that there is already some loss of correlation when the comparison begins. This effect 
merely means that the graph for the exponential function for these sample lines 200, 210, 220, 
230 is shifted to the left by an amount proportional to the initial value for loss of correlation; it 

30 has no impact on the computation of the angles of separation 0\ or 0 2 . 
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[0036] Similarly, the position of planes which are parallel to each other in either or both of the 
(0,0,0) to (p, 0, 0) and (0,0,0) to (p, 0, tt/2) directions can easily be computed using the methods 
discussed above. For any direction in which a plane is parallel to the reference plane 100a, the 
rate of correlation loss will be zero, and thus the angle of separation will be zero. If the plane 
5 does not intersect the reference plane 100a, then the initial value of the correlation loss function 
will be non-zero, but rate of correlation loss will remain at that non-zero value, thus indicating a 
parallel but non intersecting plane, in the direction of the correlation loss measurement. 

[0037] Thus, the position of any plane may be determined, relative to the position of any 
arbitrarily selected reference plane, whether the plane intersects the reference plane or not, 
10 whether the plane is parallel to the reference plane or not, by comparing the image data 

contained in the two planes, and computing the rate of correlation loss in each of two directions, 
or sometimes in one direction, if it is known that the two planes intersect. 

[0038] Expanding on the principles discussed above, once the position of a first plane is 
determined relative to a reference plane, then the position of a second plane relative to the first 
15 plane can be determined, by using the first plane as the reference plane and performing the 
position determination again. Thus an arbitrarily long chain of planes, positioned in three- 
dimensional space, may be constructed using the methods disclosed herein. Where these planes 
each contain image data gathered from a catheter, as the catheter travels through a lumen in a 
patient, this chain of planes represents a map of the lumen, in three dimensions. 

20 [0039] In the foregoing specification, the invention has been described with reference to specific 
embodiments thereof It will, however, be evident that various modifications and changes may 
be made thereto without departing from the broader spirit and scope of the invention. For 
example, the planes may have any orientation with respect to each other, not merely the 
orientations described above. The data compared for correlation loss could be any data for 

25 which a relative position computation is desired, and not merely the spatial image data described 
above. For example, images gathered at different points in time could be compared to determine 
age-triggered correlation losses. Further, the reader is to understand that the specific ordering 
and combination of process actions described herein is merely illustrative, and the invention can 
be performed using different or additional process actions, or a different combination or ordering 

30 of process actions. Features and processes known to those of ordinary skill in the art of medical 
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devices may similarly be incorporated as desired. Additionally, features may be added or 
subtracted as desired. The specification and drawings are, accordingly, to be regarded in an 
illustrative rather than restrictive sense, and the invention is not to be restricted or limited except 
in accordance with the following claims and their legal equivalents. 
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